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It was shown in a previous work that the data combinations canceling laser frequency noise 
constitute a module - the module of syzygies. The cancellation of laser frequency noise is crucial 
for obtaining the requisite sensitivity for LISA. In this work we show how the sensitivity of LISA 
can be optimised for a monochromatic source - a compact binary - whose direction is known, by 
using appropriate data combinations in the module. A stationary source in the barycentric frame 
' appears to move in the LISA frame and our strategy consists of coherently tracking the source by 
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appropriately switching the data combinations so that they remain optimal at all times. Assuming 
that the polarisation of the source is not known, we average the signal over the polarisations. We 
find that the best statistic is the 'network' statistic, in which case LISA can be construed of as two 
independent detectors. We compare our results with the Michelson combination, which has been 
used for obtaining the standard sensitivity curve for LISA, and with the observable obtained by 
optimally switching the three Michelson combinations. We find that for sources lying in the ecliptic 
plane the improvement in SNR increases from 34% at low frequencies to nearly 90% at around 20 
mHz. Finally we present the signal-to-noise ratios for some known binaries in our galaxy. We also 
show that, if at low frequencies SNRs of both polarisations can be measured, the inclination angle 
I of the plane of the orbit of the binary can be estimated. 
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PACS numbers: 04.80.Nn, 07.05.Kf, 95.55.Ym 



I. INTRODUCTION 

The goal of the LISA mission Q is to detect and analyze low frequency gravitational signals essentially in the range 



CT from 0.1 mHz to 0.1 Hz. In this frequency range, sources of gravitational waves are mainly the wide population of 
galactic compact binaries, and interactions between black holes with a range of complexity. The study of the emission 
j^Tj, of GW from known binaries could be extremely useful for firstly, direct determination of distances, and secondly, 
possible small general relativistic effects, if the signal to noise ratio (SNR) is large enough. For this reason, we focus 
on optimising the sensitivity of LISA for a given astrophysical source with known direction. The sensitivity of LISA 
can be improved by solving technological problems, but it can also be improved by employing certain optimal data 
analysis strategies. In this paper, we show how algebraic methods previously developed in [2J, henceforth referred 
to as paper I, can be used to design optimal strategies for combining data with appropriate time-delays. We show 
in this paper that it is possible to maintain the optimality during the year by continuously updating the parameters 
of the combination or simply by switching to optimal data combinations as the source appears to move in the LISA 
frame, as LISA moves in its complex orbit around the sun. The problem of optimisation of SNR, in this context, has 
been addressed before in Q and 3 • In Q optimisation has been carried out before averaging over the directions and 
polarisations, while in £|, henceforth referred to as paper II, the averaging is done first and then the optimisation. In 
this paper, the averaging over the polarisations is performed first and then the SNR is optimised for the average signal 
for a given direction over the relevant data combinations - those data combinations canceling laser frequency noise. 
Thus in this optimisation, the direction of source is assumed to be known, but not its polarisation. This would be 
the case for several binaries in our galaxy. We analyse two observables in this context: (i) optimal data combination 
in the module which yields the maximum SNR for a given direction and (ii) a 'network' observable which is obtained 
by squaring and adding the SNRs of two independent (orthogonal) data combinations one of them being the optimal 
combination mentioned in (i) and another orthogonal to it. The network observable in general yields higher SNR. We 
analyse how these SNRs depend on direction of the source. For an integration time of one year, we compare our results 
by plotting sensitivity curves with those obtained from, (a) the Michelson combination X which has been used 
for obtaining the standard sensitivity curve for LISA, and (b) the observable A" sw i tc h obtained by switching the three 
Michelson combinations optimally. We find that for sources lying in the ecliptic plane, the network sensitivity at low 
frequencies is about 34% more than the optimally switched Michelson combinations which rises to nearly 90% at 20 
mHz. Finally we compute the optimal SNRs for six known binary systems in our galaxy, whose SNRs are significantly 
high, for an integration time of one year. We show that if at low frequencies the two SNRs of the orthogonal data 
combinations can be measured, then it is possible to estimate the inclination angle of the binary's orbit. 
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FIG. 1: The schematic LISA configuration 



-v 




The paper is organised as follows: In section [H] we briefly describe the formalism involving commutative algebra 
developed in paper I. We also present a convenient set of generators for the data combinations found in paper II, which 
consisted of the eigenvectors of the noise-covariance matrix, and in which the computations are simplified. In section 
IIIII we discuss the gravitational wave (GW) signal from non-spinning compact binaries and average the signal over 
the polarisations. In section Hvl we obtain the optimum SNR by the method of Lagrange multipliers, where we must 
solve an eigenvalue problem. In section we describe the transformations from the LISA frame to the barycentric 
frame and vice- versa. In section IVT1 we compute the optimal SNRs and the sensitivity curves and compare them with 
the standard curve for LISA obtained with the Michelson data combination X and the optimally switched Michelson 
combination X sw j tc h. In this section we also estimate the SNRs for a few known binaries in our galaxy and describe 
a method to estimate the angle of inclination of the binary's orbit, if the two SNRs of the orthogonal combinations 
can be measured. 



II. THE MODULE OF SYZYGIES AND ITS GENERATORS 

In this section, we briefly review the results from papers I and II which arc needed for the analysis that follows. In 
LISA the six elementary data streams are labeled as U l and V 1 , f = 1,2, 3; if the space-crafts are labeled clockwise, the 
U l beams also travel clockwise, while the V 1 beams travel counter-clockwise. The beam U 1 travels from space-craft 
3 to space-craft 1 along the arm of length L2 in the direction — A2, while — V 1 represents the beam traveling from 
space-craft 2 to space-craft 1 along the arm of length L3 in the direction of 113. The remaining 4 beams are described 
by cyclically permuting the indices. These beams contain the laser frequency noise, other noises such as optical path, 
acceleration etc. and also the GW signal. For performing the analysis, we choose a frame (xl,Ul, z£) tied to the 
LISA constellation. We choose the centroid of the LISA triangle as the origin, the LISA triangle to lie in the (xl,Ul) 
plane and space-craft 2 to lie on the a^-axis. We will henceforth drop the suffix l L' from the LISA frame quantities 
and only introduce it when other frames are being considered in the discussion. A schematic diagram of LISA is 
shown in FigQ] 

In paper I, we showed that all data combinations canceling laser frequency noise form a module over the ring of 
polynomials in three indeterminates, these being the three-time delay operators Ei, i = 1, 2, 3 of the light travel time 
along the three arms with lengths Li respectively. Thus for any data stream X(t) : EiX(t) — X(t — Li). For the 
LISA specifications Li ~ 16.7 sec. corresponding to an armlength of 5 million km. A general data combination is a 
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linear combination of the elementary data streams U l ,V l : 

3 

X(t) = ^2p i V i (t)+q i U i (t) ) (1) 

i=l 

where p t and q,; are polynomials in the time-delay operators E^i = 1,2,3. Thus any data combination can be 
expressed as a six-tuple polynomial 'vector' (j>i,qi). For cancellation of laser frequency noise only the polynomial 
vectors satisfying this constraint are allowed and they form the module of syzygies mentioned above. An important 
advantage of this formalism is that, given the generators of the module, any other data combination canceling the 
laser frequency noise is simply a linear combination of the generators with the coefficients being polynomials in the 
ring. 

Several sets of generators have been listed in paper I. Based on the physical application as also on convenience, we 
may choose one set over another. For the monochromatic sources, as argued in paper II, only three generators are 
sufficient to generate the relevant data combinations. In general 4 generators are necessary to generate the module 
but if the source is monochromatic, except at certain frequencies which are solutions of e l ( L i+ L 2+ L z) n = i ) the fourth 
generator can be effectively eliminated for the purposes of maximising the SNR. Since this maximisation is possible 
arbitrarily close to the singular frequencies the singularities do not seem to be important. 

We found in paper II that a convenient set of generators for our purpose is the one that diagonalises the noise 
covariance matrix. These generators are the eigenvectors of the noise covariance matrix. The noise covariance matrix 
Af^jyl, J = 1,2,3 is defined as the outer product ■A/j-jj = N^Nfo where is a 12 dimensional complex noise 

vector of the generator = (yPi\qf^ (The could be any set of generators not necessarily the eigenvectors 

of the noise covariance matrix J^f^jl )• The noise vector is given by: 

= (y^(2pf ] +rl I) ),V^(2q[ I) +4 I) ),V^p[ I \V^ql I) y (2) 

where the polynomials are defined through the equations r% = — (pi +£'392) = — (<Zi +-E2P3) plus cyclic permutations 
for r 2 and r 3 . The S pf = 2.5 x 10- 48 (//lHz)- 2 Hz~ 1 and S opt = 1.8 x 10- 37 (//lHz) 2 Hz~* are the one-sided power 
spectral densities (psd) of the proof-mass noise and the optical-path noise respectively P, 0] . 

For the purposes of this paper when computing the SNR for any given data combination, the differences in arm- 
lengths are extremely small compared to the wavelength of the GW we are interested in detecting. Thus we may 
take the armlengths to be equal, i.e., E\ — E 2 — E 3 — E = e lQjL . For the purposes of this analysis it is sufficient 
to consider LISA as an equilateral triangle; we ignore the deviations arising from this assumption. This simplifies 
the expressions for the noise and further also for the signal. The eigenvectors now contain overall common factors 
which are polynomials in E. These common factors make no difference to the computation of SNR as they cancel 
out from the numerator and denominator which comprise of the signal and the noise respectively. The unnormalised 
eigenvectors (with common factors canceled) of the noise covariance matrix is the set {l^ 7 -'} which we list below: 

Y" (1) = (1 - E, 1 + 2E, -2 - E, 1 + 2E, 1 - E, -2 - E) , 

y (2) = {-E - 1, 1, E, -1, 1 + E, -E) , (3) 

y (3) = (1, 1, 1, 1, 1, 1). 

In paper II, the have been listed with the uncanceled common factors. 

We adopt the following terminology: we refer to a single element of the module as a data combination; while a 
function of the elements of the module, such as taking the maximum over several data combinations in the module 
or squaring and adding data combinations belonging to the module, is called as an observable. The important point 
to note is that the laser frequency noise is also suppressed for the observable although it may not be an element of 
the module. 



III. THE GW SIGNAL FROM BINARIES 
A. The waveform 

Since binaries will be important sources for LISA the analysis of such sources is relevant. One such class is of 
massive or super massive binaries whose individual masses could range from 10 3 M Q to 10 8 Mq and which could be 
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around a Gpc away. Another class of interest are known binaries within our own galaxy whose individual masses 
are of the order of a solar mass but are just at a distance of a few kpc or less. Here we will assume the direction of 
the source to be known, which is justified for known binaries in our galaxy; but even for the former case of distant 
binaries, it amounts to 'looking' in a specific direction. 

The spacetime metric perturbation for a gravitational wave propagating in a w direction can be written as, 



hij (t, r) = h+(t-w-r) (9 l 9 j - fafc) 

+h x (t - w ■ f) [Oifa + 6j<t>i) , 



(4) 



where the direction to the source is described by direction angles (0, 0), the unit vector pointing to the source is 
w = (sin 9 cos <f>, sin 9 sin cf>, cos 9) and 9 and <f> are unit vectors in the direction of the coordinates 6 and <f>, such that 
{u>, 9, <fi} form a right-handed triad. 

For a binary that is not chirping, and hence monochromatic with frequency f2, the two polarisation amplitudes 
h + (t) and h x (t) Q are given by 



h+(t) = H 
h x (t) = H 



1 + cos e 



2 

1 + cos 2 e . 



cos 2ip cos fit + cos e sin 2tp sin fit 
sin 2i[> cos fit + cos e cos 2ip sin fit 



(5) 



where the angles (e, ip) describe the orientation of the binary orbit (e, ip could be the direction angles of the orbital 
angular momentum vector) and H the overall amplitude which depends on the masses, the distance and the frequency 
/ as given below: 



H = 1.188 x 1(T 22 
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(6) 



where R is the distance to the source, / = Vl/2tt is the GW frequency of the source and M. = (/i 3 M 2 ) 5 the chirp-mass, 
where [i and M are respectively the reduced and the total mass of the binary. For the typical parameters taken above 
the frequency evolves very slowly, so much so that the time for the system to coalesce is more than 20 years. For an 
observational time of the order of an year, the frequency of the binary changes very little so that the source can be 
practically taken to be monochromatic. 

Since we deal with essentially monochromatic sources, the Fourier domain is appropriate for further analysis. In 
the Fourier domain we have, 



h+(Q) = H 
h x (Q) = H 



1 + cos e 
2 

1 + cos 2 e 



cos 2ip — i cos e sin 2ip 
sin 2ip — i cos e cos 2ip 



(7) 



B. The signal matrix 

The GW response for a generator = (pf\ q^P) is (see papers I and II), 



where, 



(8) 



Fi I} (n) = ]T [p^Fv^ + q\ I] F m ., + ) (0) , 

i=l 

i#>(n) = ^(p^Fy^+q^Fu^ifl), 



(9) 
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where we have expressed the response of the two polarisations in terms of the responses of the elementary data streams 
for each of the polarisations. Below we state the responses for the + and x polarisations for the beams U 1 and V 1 
only; the responses for the remaining four beams U 2 , V 2 , U 3 , V 3 are obtained from cyclic permutations: 

ifl(w-f 3 +L 2 ) , . 

Fu M = ^JT—— T (l-e- <ni2(1+ *' ft2) W,x, 



where, 



2(1 + w ■ h 2 ) 

in(w-f 2 +L 3 ) , , 

Vx;+,x = ~ ;— ^ (l-e-^3(l— n 10 

2 (1 — w ■ n-i) V / 



6;+ = {6 -hi) 2 -(</>■ hi) 2 , & ;X = 2(0-hi){<f>-hi). (11) 
We define the signal covariance matrix h/1 in analogous fashion as the noise covariance matrix. It is given as follows: 

h (I) - h^h* 
(./) - n (J) ' 

= (>j J) /i + + if ) /i x ) (F +{ j ) h + + F x{J) h x y . (12) 

In general we may not have any knowledge of the polarisation of the GW binary source. We therefore average over the 
polarizations and assume that the direction of the orbital angular momentum of the binary is uniformly distributed 
over the sphere. The orientation of the binary (its orbital angular momentum vector) has been described in terms of 

the angles e and ip in Eq.Q. Thus we carry out the averaging of over (e, tji) which results in an overall factor of 



V) 

2/5. The averaged matrix we denote by In the Fourier domain it is given by, 



W$(fi) = H 2 (0 (if ^ + F^F* (J) ) (n) . (13) 

The signal matrix so averaged has the following properties: 

• TL is the sum of outer products of two vectors: with its complex conjugate and F x ^ with its complex 

conjugate. Thus the natural basis for expressing Ti consists of the two vectors and F x \ In the analysis 
that follows we will use this fact. 

• Because we average over the polarisations, TL is constructed out of two vectors. Its rank is two, everywhere 
except on the 9 = § plane where it is one when F$P goes identically to zero. In paper II we had obtained 
a signal matrix of rank three because there we had averaged over the directions as well. While in since 
optimisation is performed first before averaging, the signal matrix is constituted from a single vector and thus 
has rank one. 



IV. OPTIMIZING SNR 

For a generic data combination a^Y^ where am are polynomials in E, the SNR is given by: 



SNR 2 = (14) 



If the psd is given in units of Hz -1 then the above equation yields the square of the SNR integrated over one second. 
Because one second is short compared to the various time-scales envisaged in the problem we call this SNR the 
instantaneous SNR. Since we are dealing with monochromatic sources, E = e , the coefficients am reduce to just 
complex numbers. Also in the generating set {^ < ' / ' ) } the noise covariance matrix is diagonal with diagonal elements 
iVjy However, we find nm = n^ 2 ) m our case, so that the first two eigenvectors correspond to the same eigenvalue. 

Thus, the denominator of Eq. l|14fl simplifies to a sum of squares \ct^x) | 2 n i + \ a (2) | 2 n i + \ a (3) | 2 n \ which we can set 
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equal to unity because the SNR does not depend on the normalisation of the data combination. Moreover, it is 
convenient to define coefficients which are scaled by the noise, 0m — ctmn\, 0i 2 ) = a (2) n 2, 0(3) — "(3)^3 so that 
the 0m satisfy, 

|/?(l)| 2 + |/?(2)f + |/? (3 )| 2 = 1. (15) 

The expression for SNR simplifies to, 

SNR 2 = a T -H-a = T -p-0, (16) 

where we define a SNR matrix p by, 

P ^ ~ n (I) n (J) ■ (17) 

In the above equations the a, are construed of as 3 x 1 column matrices and H and p as 3 x 3 square matrices. 

The extremisation of SNR is now carried out with the method of Lagrange multipliers because of the normalisation 
constraint on 0. This procedure yields an eigenvalue equation with the Lagrange multiplier appearing as an eigenvalue: 

p-0 = X0. (18) 

Since TC has atmost rank 2, one eigenvalue is necessarily zero. We now proceed to compute the other two eigenvalues. 
The analysis is simplified if we go to the basis consisting of the two vectors /+ and f x ■ 

tf ] =h ^, fi T) =h ^, (19) 
nm n {1) 

where h Q = (yj2/h)H is the amplitude of the GW averaged over the polarisation states at frequency VI. We can then 
write the matrix p as a tensor product in terms of these two vectors: 

p = f + ®f++fx ®/x*- (20) 

In general the vectors /+ and / x are not orthogonal. 
The action of the matrix p on any vector v is given by, 

p.v = (/+ ® /; + /x ® /* ) • v = (/; • v)f + + (/* • v)f x . (21) 

For the eigenvalue problem, we have the eigenvalue equation, 

(/; • v)f+ + (f* ■ v)f x - Xv. (22) 

Expressing the eigenvector v as a linear combination of / + and / x , v = c + / + +c x fx , and from the linear independence 
of /+ and f x (they are linearly independent in general) we obtain the system of equations for c + and c x as: 



(|./+| 2 -A) c+ + (/;./ x ) Cx = 

(/+-/x)c + + (|/x| 2 -A)c x = 0. 
Setting the determinant of this system to zero, we obtain eigenvalue equation: 



(23) 



A 2 -(|/+| 2 + |/x| 2 )A+|/ + x/ x | 2 = 0, (24) 
which can be solved to yield the eigenvalues: 

A± = i(|.f;| 2 + |/x| 2 ±A/A), (25) 

where, 

a = (i/;i 2 + i/ x i 2 ) 2 -4|/; x / x | 2 
= (i/;i 2 -i/xi 2 ) 2 +4i/;-/ x ?. (26) 
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In the low frequency limit as we shall find /+ • / x ~ so that the positive square root of A is just |/+| 2 — |/ x | 2 and 
the eigenvalues are just |/+| 2 and |/ x | 2 (infact the eigenvalue equation factorises into linear factors) with / + and / x 
as eigenvectors respectively. This can also be directly inferred from the structure of p. So it may be appropriate to 
call A_ as A x even in the general case. In the general case the eigenvectors can be easily determined from Ea. (|23|) : 

v+ = (a+-i/;i 2 )/; + (fV-f+)fx, 

v, = (fx ■ /+*)/+ + (A x - \f+\ 2 )fl , (27) 

where u+. x are eigenvectors belonging to the eigenvalues A+ iX respectively. Since p is Hermitian the eigenvectors are 
orthogonal, v + ■ v* x = as can be verified also directly from the expressions of the eigenvectors. These eigenvectors 
are not normalised. 

The eigenvalues A+, x are the squares of the instantaneous SNRs for the two data combinations described by the two 
corresponding eigenvectors. The data combinations are v + n\Y^ and v y ^Y^, which we will call eigen-combinations 
or alternatively eigen-observables, and which give the instantaneous SNRs: SNR 2 ^ x = A +iX respectively. For a given 
direction (9, <j>) in the LISA frame, SNR + is the maximum instantaneous SNR among all data combinations. While 
SNR X is the minimum instantaneous SNR among data combinations that are linear combinations of v + and v x ; those 
which lie in the 'plane' of v+ and ? x . However, SNR X is not zero in general and therefore not the absolute minimum; 
the absolute minimum SNR is zero corresponding to the third eigenvalue which is zero. Moreover, the eigenvectors are 
orthogonal, which means they yield statistically independent observables. Thus these observables can be combined in 
quadratures to form a network observable with instantaneous SNR ne t WO rk given by: 

SNR-nctwork = SNR+ + SNR 2 , 

= A+ + A x =|/;| 2 + |/;| 2 . (28) 

The third eigenvector is /+ x f x ; it is orthogonal to / + and / x with eigenvalue zero Qjj . This means that the data 
combination corresponding to this vector gives zero response in that particular direction, which may be important if 
one wishes to 'switch off' the GW coming from that direction. 

V. TRACKING A GW SOURCE WITH LISA 

In general the amplitude of the GW source will be small and it would be necessary to track or follow the source 
for a considerable period of time in order to accumulate an adequate SNR. This period of time could range from few 
days to a year or even years - the life of the LISA experiment. The LISA configuration performs a complex motion 
due to which the source will appear to move in the LISA frame even if it is stationary in the barycentric frame. For 
the purposes of this analysis, we will take the observation time to be a year and integrate the SNR for this period 
of time. It is possible that some GW sources may be sufficiently powerful that SNR integration for an year is not 
required. In that case it is easily possible to limit the integration to the required period of observation and obtain 
useful results. We will also assume that the source is stationary in the barycentric frame, that is, its direction remains 
constant during the observation period. In this section we will present the transformations connecting the barycentric 
frame and the LISA frame and hence obtain the apparent motion in the LISA frame for a GW source fixed in the 
barycentric frame. 

The LISA constellation trails the Earth in its orbit by 20° around the sun. The plane of LISA makes an angle 
of 60° with the plane of the ecliptic and the LISA triangle rotates in its own plane completing one rotation in a 
year. We describe the barycentric frame by the Cartesian coordinates {xb,Vb, zb}- The xb — Vb plane coincides 
with the orbital plane of LISA. The ze-axis is orthogonal to this plane forming a right-handed coordinate system. 
The transformation from the barycentric frame to LISA frame is given as follows 0: A vector re = (xb, Vb, zb) T is 
transformed to tl = (xl, Vl> zl) T by the matrix TZ as, 

r L =ft-r B , (29) 
where 1Z is a product of the three matrices containing Euler angles: 



1l = C-B- A, 



(30) 
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where the matrices A, B and C are given by, 



cos ip a sin ip a 
A = \ — sin ip a cos ip a 
1 



B 




cos ipc sin ip c 

— sinip c cosV'c | . (31) 
1 

Here ipa — ojt + ao, ip c = —ojt + flo and u> — 2tt/Tq. Tq is the orbital period of LISA which we take to be of one 
year duration. The ao and /3 are constants fixing initial conditions when the observation begins. The matrix 1Z is 
time-dependent and is a product of two time-dependent matrices. 

The unit vector wb of the source direction described by the angles (9b,<Pb) in the barycentric frame is given by 
wb = (sin#e cos </>b , sin 6*b sin 4>b, cos 9b)- The corresponding vector in the LISA frame is then, = lZ(t) ■ wb- 
From this vector we can explicitly work out the angles (9l, 4>l)'- 




(32) 



where, 



w\ = cos ip c sm. 9 B cos (4> B - ip a ) 

+ i sin ip c sin 9 B sin (cp B - tp a ) 

j.^ • i a 
+ — sm %p c cos 9 B , 

w\ = - sin ip c sin 9 B cos (cpB - ipa) 
+ icos7/> c sin6» B sin (<?!>£> - ip a ) 

V3 

+ -^- cos ip c cos 9b- (33) 

For a fixed source direction (9b,<Pb) in the barycentric frame, the apparent source direction depends on time in 
the LISA frame; 9l,4>l ar e functions of time; the source appears to move in the LISA frame. Thus a given data 
combination even if it is optimal initially, will not continue to remain optimal subsequently. Our strategy is then to 
switch the data combinations continuously so that the SNR remains optimal at all times. In this way we can optimally 
track the source and accumulate maximum SNR. The following figures show the apparent trajectory of the source 
for one year in the LISA frame when (i) the source lies at the pole of the barycentric frame - it is just a circle in the 
LISA frame: 9l = 7r/3 (See Fig|2J) and (ii) the source lies in plane of the LISA orbit - a figure of 8 is described by the 
source (See Fig|3J). 



VI. OPTIMAL SENSITIVITIES 



For optimal tracking of the source, we switch the data combinations so that the SNRs yielded are optimal at 
all times. Three SNRs are of interest; SNR +yX and SNR nc t wor k, that is, the two eigenvalues and their sum. The 
eigenvalues are functions of 9l,4>l which are in turn functions of t as the GW source describes its apparent motion 
in the LISA frame. The three integrated SNRs are given by: 

SNR^ x (w B ) = / \ +>x (6 L {t;w B ),<p L {t;w B ))dt, 
Jo 

SNR^ etwork (^ s ) = SNR^(w B ) + SNR 2 x (^ B ), (34) 
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FIG. 2: Apparent position of the source in the sky as seen from LISA frame for (6b = 0, <f)B =0). The track of the source for 
a period of one year is shown on the unit sphere in the LISA frame. 




FIG. 3: Apparent position of the source in the sky as seen from LISA frame for (8b = f , 4>b = ). The track of the source for 
a period of one year is shown on the unit sphere in the LISA frame. 
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where we have taken the total observation time to be T and the initial time of observation to be zero. If we integrate 
for a complete year, we can set the initial observation time to be zero without loss of generality and at the same time 
set the constants ao, [3q to be zero. However, if we make observations for times that are not integral number of years, 
the constants cvrj;A) must be chosen appropriately in the transformation matrix 1Z. 

An important point is that, since LISA moves in an heliocentric orbit, a Doppler phase depending on the position 
vector of the LISA projected on to the direction to the source wb and the GW frequency will be added to the phase 
of the GW signal. This Doppler phase will be a function of time which will be added to the monochromatic part of 
the phase Qt of the signal. We assume in this analysis that this phase has been accounted for when integrating the 
signal. One way is to 'stretch' the time-coordinate so that the signal appears monochromatic (a technique well-known 
to radio astronomers) [t| 0] . 

A. The Low Frequency Limit 

An important case arises when we consider GW of low frequencies, say below 3 mHz. For this case it is possible to 
obtain analytical expressions for the optimal SNRs. A large fraction of the sources for LISA fall into this category, 
for example, massive/supermassive blackhole binaries, several galactic and extragalactic binaries which contribute to 
the 'confusion noise'. 

We need to compute the /+ and / x in order to compute the SNR matrix p and then obtain its eigenvalues. 
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Integrating the eigenvalues for the period of observation will yield the relevant SNRs. To this end, we expand Fui-,+, x 
and Fvi;+,x to the lowest order in the dimensionless frequency f2L: 

FlJl;+,X = i-^ (1 + i^LT 2 )^2;+,x, 

Fvi;+,x = -i^(l + iflLT 3 )^ ;+x , (35) 



1 



2 V \/3 

^l+COS 2 6>) / 7T\ 1. 2 



where, 



£m;+ = ^ cos {2(f)- (2m- 1)-J - -sin 

Cm ; x = -cosflsin (2(f>- (2m- 1)^ , (36) 

where, f m is the unit vector in the direction of m-th spacecraft and m — 1,2,3 and the angles 9, <j) refer to the LISA 
frame (as before we have dropped the subscript 'L'). The transfer functions for the four other elementary data streams 
are obtained by cyclic permutations. In order to get the /+, / x we must operate on the U l , V 1 with the polynomial 

operators pf\ given in Eq.Q) and then scale them by the averaged signal divided by the noise. Thus, we obtain 
after some algebra, 

(l + cos 2 0), 
f+ = Pa- ^ -(-bih^cobS.Q), 

fx = -pocos^(cos$,sin$,0), (37) 

where, 

Po = ±*(nL) 2 , * = 2^+£. (38) 
V5 n\ 3 

The vectors / + and / x are expressed in the basis and in this basis they have real components. We observe the 
following properties of the vectors: 

• The vectors /+ and f x he in the (Y^',Y^ 2 ') plane and the Y^ component is zero for both vectors. This can 
be understood if we recall that Y^ is proportional to symmetric Sagnac combination Q which is insensitive to 
GW at low frequency. 

• The apparent motion of a GW source in the LISA frame can be optimally and continuously tracked by /+ and 
f x by rotating the pair (/ + , / x ) by <I> at each instant of time as follows: 



fx \ _( cos $ sin $ \ / Y (1) 
/ + J _ ^ - sin $ cos $ J y y ( 2 ) 



(39) 



Here /+ and f x are unit vectors in the directions of /+ and f x respectively. Thus, optimally tracking a source 
amounts to orienting the data combinations along / + and f x ■ 

• Moreover, the vectors are orthogonal: f + ■ f x = in this low frequency limit. The orthogonality implies that 
these eigen-observables /+ and f x have zero response to x and + polarizations of GW respectively. We use this 
fact to estimate the polarisation angles; namely (e, tp) at the end of this section. 

The eigenvalues of the eigen- vectors /+ and / x are just |/+| 2 and |/ x | 2 and are explicitly given by, 

, ,'l + cos 2 6l x 2 



A+ = Po , 2 

A x = p 2 cos 2 6i. (40) 

The eigenvalues are the squares of the instantaneous SNRs. We notice that po is the maximum instantaneous SNR 
obtained when the source lies at the poles 9 = or ir in the LISA frame. If the source is observed over a period of an 
year, the eigenvalues must be integrated over this length of time. We notice that the eigenvalues do not depend on <fi. 
Thus our next task of integrating the SNR becomes somewhat simplified. 
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FIG. 4: Instantaneous SNR+ tX and SNR nc t W ork as functions of time for the source direction (9 b = f , <f>B = 0) at the GW 
frequency of / — 1 mHz in units of po- 
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1. The integrated SNR 

For a given source direction (9b,4>b) the corresponding track (^z,(t), <j>L,{t)) of the source in LISA frame is given 
by Eq.(J22J). We may substitute these values into the integrals for the SNRs and the integrals yield simple analytical 
expressions for the SNRs if the integration is over an integral number of years. The instantaneous SNRs however 
change with time. In the FigQJwe show how the SNRs change with time as LISA orbits the sun during the course of a 
year. The various SNRs are shown for a source lying in the ecliptic plane, (9b = "ft/2, <f>B = 0) for the GW frequency 
of 1 mHz. We have chosen this direction because the SNRs show considerable variations during the course of a year. 

Integration over a period of T — Tq leads to the following results: 

SNR + (0 b ,<I>b) = SNR g+(cos9 B ), 

SNR x (9 B ,<t> B ) = SNR 5x (cos0 B ), (41) 

where, 

snRo = poV^=4f— mfVr^, (42) 

V5 nx 

and 




We have purposely not 'simplified' the formulae in powers of x 2 because in this form it is easy to see the limits 
x = ±1,0 corresponding to 9b = 0,7r,7r/2 respectively. Infact since only x 2 occurs in the expressions of g+,x there is 
symmetry about the ecliptic plane. The network SNR is just the root mean square of the two SNRs: 



SNR nctw ork(#B, 4>b) — SNR ^network 

(cos#b) , 



(44) 
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FIG. 5: The functions g+,x and (/network as functions of 9b 
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where, ^network! 2 -) = 9+( x ) +5x( a; )- I n FiglSJwe plot the functions <7+, x and (^network as functions of Ob between 
0° < 0b < 180°. The factors g are of the order of unity and the SNR gives essentially the integrated SNR of a GW 
source. Recall that this is an SNR averaged over the polarisations. 

We observe from FigJSJthat the maximum integrated SNR is obtained for sources lying in the ecliptic plane {0b = 
90°). This can be readily explained from Fig|3| where the trajectory of such a source is plotted. One observes that 
a large fraction of the orbit in the LISA frame is away from the LISA plane (0l — 90°). As seen from Ea . 1)40(1 the 
sensitivity of LISA increases as we go away from the plane of LISA, that is, towards the poles. We also observe that 
the variations in the three curves are small; among these the network SNR shows the highest variation of ~ 15%. 
This shows that LISA has a more or less uniform average response over the year as it moves in its orbit. The above 
formulae are also more generally valid if the integration time T is taken to be in integral multiples of Tq ; then the 
T Q in Eq.O should be replaced by T. 

We compute the SNRs of six binaries in our galaxy which give high average SNR (averaged over polarisations) if 
we integrate the SNR optimally over the period of a year. The following table lists six binary systems which give 
the network SNRs ranging from about 3 to over 100. The SNRs have been computed assuming circular orbits for 
the binaries. The information about the binaries has been obtained from 8]. We observe that the binary masses are 
small; ~ .5M® and companion mass ~ .02M®. The reason for such small masses is that the orbital periods of the 
binaries must be short, in this case ranging from about .01 to .03 of a day. The GW quadrupole frequency f gw is 
related to the orbital period by the equation: 

Thus the sources radiate GW at frequencies ~ 1 or 2 mHz, the band in which LISA has maximum sensitivity. This 
yields the high SNRs. They are also close by; R ~ 100 or 200 pc which tends to increase the raw gravitational wave 
amplitude H . Note that the SNRs listed in the table are average and actual observations can yield different values 
depending on the orientation of the binary orbit. Also the observables used are optimal in the average sense. If the 
orientation of the binary is known then in general a better observable can be found. 



2. Estimation of inclination angle of the orbital plane 

We showed in previous sections that if we do not have any prior information of the inclination angle e of the binary 
orbit and the the angle ip, then tracking the source with /+ and f x are optimal in the average sense. Typically, the 
orbital inclination is difficult to estimate from other astrophysical observational means. However, for binaries with 
known masses and distances (e.g. binaries in table-QJ, we can estimate (e, tp) from the output of eigen-observables /+ 
and f x . 
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TABLE I: In this table we list six binaries in our galaxy whose parameters are known. For these we compute the optimum 
SNRs where the optimisation has been performed after averaging over the polarisations. 



Name 


P or b in days 


mi 


in? 


log(tf) 


e° B 


R in pc 


SNR+ 


SNR X 


SN Rnetwork 


AM CVn 


0.011907 


0.5 


0.03 


-21.4 


124.46 


100 


89.9 


75.6 


117.4 


CP Eri 


0.019950 


0.6 


0.02 


-22.0 


116.43 


200 


8.9 


7.4 


11.4 


CR Boo 


0.017029 


0.6 


0.02 


-21.7 


72.103 


100 


26.3 


22.7 


34.7 


GP Com 


0.032310 


0.5 


0.02 


-22.2 


66.997 


200 


2.2 


1.8 


2.8 


HP Lib 


0.012950 


0.6 


0.03 


-21.4 


85.040 


100 


77.6 


67.8 


103.0 


V803 Cen 


0.018650 


0.6 


0.02 


-21.7 


120.32 


100 


20.6 


17.5 


27.0 



The integrated SNRs over the period of one year of /+ and / x , without averaging over polarisations are given by 

SNR+(0 B ,<l> B ) = y|sNR <;+ (cos B ) a+ , 

SNR x {9 b ,4>b) = y|sNR 5x (cos# B ) a x , (46) 

where a+ = \h+(Q) / H\ and a x = \h x (Q)/ H\ (see Eq. (J2J). We note that here the factor y/5/2 appears, since there 
is no averaging over the polarisations. The a +yX can be estimated, if the SNRs appearing on the LHS of Eq. I|46|) can 
be measured. Further, straightforward algebra shows that 

2 , 2 _ /l + cos 2 e\ 2 2 



a% + aj = ^ J + cos^ e. (47) 

From the above equation, we can estimate e. Substituting back in a + , one can also estimate ip if needed. This exercise 
can be carried out for the binaries listed in the tabled 



B. The general case 



In this section we relax the condition of dealing only with low frequencies and consider the entire band-width of 
LISA. We compare the sensitivities for LISA obtained by using the optimum SNR with the Michelson data combination 
usually denoted by A. As a polynomial vector (p^, q{) it is given by: 

X = {l-E 2 2l Q, E 2 (E 2 ElE 3 (E 2 2 ~l),0). (48) 

When we set all the Ei to be equal, the factor (1 — E 2 ) factors out and one is left with a simple polynomial 
vector (1, 0, — E, 1, — E, 0). This combination has been used to plot the standard LISA sensitivity curve. Two other 
Michelson observables are obtained by cyclic permutations called Y and Z in the literature [|J. In this section 
we will compare the sensitivities obtained by using the eigen-combination v + and the network observable with the 
Michelson combination X and the observable A sw i tc h obtained by 'switching' the Michelson's X, Y, Z optimally, that 
is, SNRx svitch = m&x(SNRx, SNRy, SNRz). Extending the definition of sensitivity of a data combination from 
earlier literature to an observable, we define the sensitivity of an observable W as, 

Sensitivity w (/ ) = ^- M , (49) 

where SNRyy(f) is the integrated SNR over a observation period T. The number 5 has been chosen following earlier 
literature. Eq. I|49|l for a fixed data combination Wq reduces to the standard one in literature: 



Sensitivity^ (/) = 5 V ' , (50) 

where B = 1/T. As before we take T = T Q for plotting the sensitivity curves in the figures below. The results of 
our findings are displayed in the plots. FigEl displays the sensitivity curves for the observables (a) Michelson - X 
(dotted curve), (b) Switched Michelson A switc h (dash-dotted curve), (c) Eigen-combination v + (dashed curve) and 
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FIG. 6: Sensitivity curves for the observables: Michelson, Switched Michelson, v+ and network for the source direction 
6 B = 90°, 4>b = 0°. 
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FIG. 7: Ratios of the sensitivities of the observables network, «+, x with X sw ; tc h for the source direction Ob = 90°, 4>b = 0°. 
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(d) network observable (solid curve). We observe that the sensitivity over the band- width of LISA increases as we 
go from (a) to (d). We observe that the X sv/ i t ch does not do much better than X. This is because for the source 
direction chosen X is reasonably well oriented and there is not much switching to Y and Z combinations. However, 
the network and v + observables show significant improvement in sensitivity over both X and X sw ; tc h. The sensitivity 
curves (except X) do not show much variations for other source directions and the plots are similar. The quantitative 
comparison of sensitivities is shown in Fig0 

Here we have compared the network, and the eigen-combinations v+,x with A sw itch. Defining: 



Ka(/) = 



SNR a (/) 
SNRx swltch (/) ' 



(51) 



where the subscript a stands for network or +, x and SNRx 3Witch the SNR of the observable X SW i tc h, we plot these 
ratios of sensitivities over the LISA band- width. We notice from the /^network curve that the improvement in sensitivity 
for the network observable is about 34% at low frequencies and rises to nearly 90 % at about 20 mHz, while at the 
same time the w+ combination shows improvement of 12 % at low frequencies rising to over 50 % at 20 mHz. Finally, 
FigO exhibits the sensitivities of the network observable over various source directions. Since the sensitivity of this 
observable is independent of 0s, we plot the curves for several values of 0b = 0,30,60,90 degrees. Also since the 
network observable possesses reflection symmetry about the ecliptic plane 9b = 90°, we do not need to plot the curves 
6b between 90° and 180°. The important observation from this figure is that not much variation in sensitivity is seen 
as all source directions are scanned. Thus the network observable integrated over the year has essentially uniform 
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FIG. 8: The sensitivity curves for the network observable for 8b = 0°, 30°, 60°, 90° 
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sensitivity to all source directions over the frequency range 10 4 — 10 1 Hz. 

VII. CONCLUSION 

We have shown how the SNR can be optimised for a GW source with known direction but with unknown polarisation. 
While obtaining the SNR the signal is averaged over the polarisations and then optimised. Because of this procedure 
we lose out to some extent on the SNR but on the otherhand it leads to robust results. The optimisation methods 
are algebraic in that one must solve an eigenvalue equation to determine the optimum SNRs. We separately deal 
with the low frequency case as it is of considerable astrophysical importance - a large fraction of GW sources are 
expected to be of this category. Secondly, it lends itself to simple analytical approximations which throw light on 
the results obtained. Lastly, we deal with the general case covering the full band-width of LISA. We have compared 
the sensitivities obtained with our strategy to those obtained in the standard way. We find that the improvement in 
sensitivity of the network observable over the X or A sw i tc h ranges from about 34% to nearly 90% over the bandwidth 
of LISA for a source lying the ecliptic plane. Finally we present a list of few binaries in our galaxy for which the 
optimal SNRs and the network SNRs have been computed. We also describe a method of extracting information 
about the inclination angle of the orbit of the binary if SWi? +jX can be measured. 
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